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Abstract We study how active-region-scale flux tubes rise buoyantly from the 
base of the convection zone to near the solar surface by embedding a thin flux 
tube model in a rotating spherical shell of solar-like turbulent convection. These 
toroidal flux tubes that we simulate range in magnetic field strength from 15 
kG to 100 kG at initial latitudes of 1° to 40° in both hemispheres. This article 
expands upon Weber, Fan, and Miesch (Astrophys. J., 741, 11, 2011) (Article 
1) with the inclusion of tubes with mag netic flux of 10 20 Mx and 10 21 Mx, and 
P-n \ more simulations of the previously investigated case of 10 22 Mx, sampling more 

Q_[ convective flows than the previous article, greatly improving statistics. Observed 

properties of active regions are compared to properties of the simulated emerging 
flux tubes, including: the tilt of active regions in accordance with Joy's Law as 
in Article 1, and in addition the scatter of tilt angles about the Joy's Law trend, 
the most commonly occurring tilt angle, the rotation rate of the emerging loops 
with respect to the surrounding plasma, and the nature of the magnetic field 
at the flux tube apex. We discuss how these diagnostic properties constrain the 
initial field strength of the active region flux tubes at the bottom of the solar 
convection zone, and suggest that flux tubes of initial magnetic field strengths of 
> 40 kG are good candidates for the progenitors of large (10 21 Mx to 10 22 Mx) 
solar active regions, which agrees with the results from Article 1 for flux tubes of 
10 22 Mx. With the addition of more magnetic flux values and more simulations, 
we find that for all magnetic field strengths, the emerging tubes show a positive 
Joy's Law trend, and that this trend does not show a statistically significant 
dependence on the magnetic flux. 
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1. Introduction 



The toroidal magnetic field responsible for the emergence of solar active re- 
gions is believed to be generated by a dynamo mechanism at or near the base 



of the convection zone (e.g. |Spiegel and Weiss} 119801 |van Ballcgooijcn, 1982; 



IMoreno-Insertis et all UMZl IGilmanl 120001 ICharbonneaul [2"0T0j) . Of essential 
importance in the understanding of solar dynamo theory, and indeed all of solar 
physics, since flux emergence also regulates solar variability, is addressing how 
active region flux tubes rise dynamically through a turbulent convection zone. 
Also of interest is identifying the dynamo-generated magnetic field strength at 
the base of the convection zone required to produce the observed properties of 
solar active regions. 

Valuable insights into the nature and evolution of rising magnetic loops in 
the solar convective envelope have been gained through the use of the thin flux 
tube approximation by a plethora of authors (e.g. Spruit 1981: Morcno-Inscrtis, 
rragPl IFerriz-Mas and Schusslerl [TQMl |Longcope and Klapper) 110071 IFanl [2TJ00|) . 
It is found that in order for these simulated flux tubes to exhibit tilt angles and 
latitudes of emergence that agree well with observed solar active regions, the 
toroidal magnetic field at the base of the convection zone needs to be in the 
range of »30 kG to wlOO kG (jChoudhuri and Gilmanl [10571 ISchiissler etaTX 
ITQM1 ID'Silva and Choudhuril 115911 |Caligari et aT\ [TOPS]) . These studies also 
indicate that the Coriolis force is responsible for several observed asymmetries 
that exist between the leading (in the direction of solar rotation) and following 
polarities of solar active regions, such as: the tilt angles in accordance with 
Joy's Law (jD'Silva and Choudhuril 119931 Caligari et al. I1995p . the apparent 
faster proper motion of the leading polarity of an emerging active region on 
the solar surface ( van Driel-Gesztelyi and Petrovay| [T990 Mor eno-Insertis et all 
1 19941 [Caligari et aT. I1995p . and an asymmetry where the leading polarity shows a 
more coherent morphology (jFan et aU IT993: Cal igari et al. [1995 ; Calig ari et al. 



1998). However, as these previous studies neglect the effects of turbulent con- 
vection on rising flux tubes, it is possible that convective downdrafts can pin 
portions of the flux tube down to the base of the convection zone, while helical 
upflows between the downdrafts can boost the rise of the flux tube ([Fan et ail 
20031 IWeber et a77l I2011|) . Therefore, resulting emerging loops may exhibit dif- 
ferent properties from those produced in the absence of convection and pose an 
intriguing topic of study. 

As of yet, the magnetic field strength at which the solar dynamo operates is 
not well known, nor is it directly accessible via observations. However, solar cycle 
dynamo models that incorporate the Lorentz force from large scale mean fields 
indicate that the magnetic field strength generated and amplified at the base of 
the convection zone is ~15 kG, and most likely cannot exceed 30 kG ( |Rempel[ 
2006aJ |Rempel| I2006b|) . Recent simulations of solar- like stars that rotate three 
times the current solar rate have shown that a rotating convective envelope can 
generate a dynamo that consists of opposite polarity magnetic wreaths in two 
hemispheres, which span the depth of the convection zone ([Brown et aU |20T0|) . 
When portions of these wreaths become strong enough, «35 kG or greater, a 
buoyant magnetic loop develops which then rises through the convecting fluid 



SOLA: wfm_paper_v0709_links.tex; 8 August 2012; 0:43; p. 2 



Comparing Simulations of Rising Flux Tubes with Observations of Active Regions 



in which it is embedded (jNelson et al.l 1201 1[) . While these dynamo-producing 
convection simulations are not meant to reproduce the solar dynamo directly, 
they do demonstrate that persistent toroidal fields can coexist with convection 
at moderate magnetic field strengths. In light of these studies, it is important to 
understand how toroidal flux tubes of weak to moderate field strengths, «15 kG 
to «50 kG, behave as they rise through a turbulent solar convective envelope. 
Some studies have been performed that investigate the buoyant rise of fully 
three-dimensional isolated flux tubes in a turbulent convective velocity field (e.g. 
IFan et al\ 120031 IJouve and Brunl I2009|) . However, due in part to the limited 
numerical resolution of these simulations, large values of magnetic flux must be 
used, which are greater than typical active region flux, in order to keep the tube 
from dissipating as it rises. 

Recently. I Weber et aU 120 111 (hereafter Article 1) incorporate a thin flux tube 
model in a (separately computed) rotating three-dimensional convective velocity 
field representative of the solar convective envelope, in an effort to study the 
effects of turbulent solar- like convection on the dynamic evolution of active region 
scale emerging flux tubes. Although the thin flux tube model cannot capture the 
possible fragmentation of the flux tube and its internal magnetic structure, it 
does preserve the frozen-in condition of the magnetic field. The thin flux tube 
model is also useful in that each simulation can be performed quickly on single 
processor desktops, as compared to the three-dimensional simulations, which 
require multi-million processor hours on massively parallel supercomputers. As 
such, the thin flux tube approach is a useful platform for a parameter space study, 
and provides a starting point for the investigation of the effects of convective 
flows on flux tubes of realistic magnetic field strengths with magnetic fluxes 
similar to those of active regions observed on the Sun. 

This article builds upon the work as presented in Article 1. In Article 1, it 
was found that as the magnetic field strength increases from 15 kG to 100 kG, 
the dynamic evolution of the emerging flux loops changes from being convection 
dominated to magnetic buoyancy dominated. For these flux tubes, the convective 
flow is found to reduce the rise time and the latitude of emergence through 
the anchoring of flux loop footpoints by downdrafts. The addition of solar-like 
convection also promotes tilt angles that are consistent with the observed mean 
tilt of solar active regions in part because of the mean kinetic helicity of the 
upflows in the simulated convection. In this article, we expand the study by 
increasing the number of simulations performed by about 3.5 times per magnetic 
field strength to improve the statistics of the results, sampling different time 
spans of the convective flows. Solar active regions exhibit magnetic flux in the 
range of 10 20 Mx to 10 22 Mx, which include ephemeral regions and pores to 
the largest scale sunspots (IZwaanl I1987[) . To study the dependence of flux tube 
evolution on magnetic flux, we also include cases with flux of 10 20 Mx and 10 21 
Mx in this study in addition to the 10 22 Mx cases. This in combination with the 
increased number of simulations per magnetic field strength results in a total of 
5940 flux tubes to analyze. As with Article 1, we consider flux tubes with initial 
latitudes of 1° to 40° in both the northern and southern hemispheres. Besides 
the observed active region properties described by Joy's Law, in this article 
we also consider additional observational diagnostics, including the scatter of 
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active region tilts and sunspot rotation rates to further constrain the initial field 
strength of active region flux tubes at the base of the solar convection zone. We 
also discuss the properties of the magnetic field at the apices of the emerging flux 
loops. In comparison to Article 1, our new study confirms that with convection 
all flux tubes with magnetic field strengths of 15 kG - 100 kG do exhibit a 
positive Joy's Law trend, and given the uncertainties that we obtain, we note 
no significant dependence of the flux tube tilt angle on magnetic flux. 

A description of the thin flux tube model and how it is incorporated with 
the spherical shell of solar-like convection employed in this study is outlined in 
Section [2] We highlight the results obtained from our simulations in Section O 
discussing flux tube rise times, latitudes of emergence, tilt angles, magnetic field 
properties, and the rotation rate of the simulated emerging loops. A summary 
of the results is given in Section H) 



2. Model Description 

The dynamics of a thin, isolated magnetic flux tube can be described by the 
thin flux tube approximation. These equations are derived from the ideal MHD 
equations, operating under the assumption that the flux tube radius is small 
compared to all other relevant scales of variation associated with the flux tube 
(e.g. |Spruit[ [TggS |Cheng| HMD [Longcope and Klapper[ \TM7\ \Faxj\[XMb. As in 



Article 1, we solve the following equations, which describe the evolution of each 
Lagrangian element of the one-dimensional flux tube: 
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where, r, v, B, p, p, T, which are functions of time t and arc length s mea- 
sured along the tube, denote respectively the position, velocity, magnetic field 
strength, gas density, pressure, and temperature of a Lagrangian tube segment 
with subscript e referring to the external plasma which are functions of depth 
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only, 1 = dr/ds is the unit vector tangential to the flux tube, and k = d 2 r/ds 2 is 
the tube's curvature vector, subscript _L denotes the component perpendicular 
to the flux tube, $ is the constant total flux of the tube, fj, is the mean molecular 
weight of the surrounding external plasma, g is the gravitational acceleration and 
a function of depth, f2o is the angular velocity of the reference frame co-rotating 
with the sun, with Hq set to 2.7 x 10~ 6 rad s _1 in this calculation, Cd is the drag 
coefficient set to unity, 7 is the ratio of specific heats, R is the ideal gas constant, 
and v e (r, t) is a time-dependent velocity field relative to the rotating frame of 
reference that impacts the dynamics of the thin flux tube through the drag force 
term. In the above equations, we do not introduce an explicit magnetic diffusion 
or kinematic viscosity term. However, we do introduce a drag force (last term 
in Equation ([T])). which describes the interaction of the external fluid with the 
flux tube in a high Reynolds number regime (e.g. [BatchelorJ I1967p . The flux 
tube evolves passively in the external fluid such that the tube imparts no back 
reaction on the fluid in which it is embedded. The magnetic field of the flux 
tube is untwisted such that it only has a component in the I direction, and the 
tube is discretized with 800 uniformly spaced grid points along its arc length 
s. A description of the numerical methods used to solve the flux tube evolution 
as determined by the above set of equations is discussed in detail by IFan et al\ 
(T19931) . 

The thin flux tube approximation is no longer satisfied when the radius of 
the flux tube is on the order of the local pressure scale height. With this in 
mind, the thin flux tube approximation is satisfied for the bulk of the convection 
zone, where the flux tube spends the majority of its time during evolution. For 
large flux values of «10 22 Mx, and small magnetic field strength of wl5 kG, the 
radius of the flux tube will be the greatest, and will expand significantly as it 
rises toward the top of the convection zone. Due to these effects, we must stop 
our simulation before the flux tube can reach the photosphere. In this case, we 
stop our simulations once the flux tube has reached a height of «21 Mm below 
the solar surface. 

A reference one-dimensional solar structure model 



( Christensen- Dalsgaar d et al. I1996P provides the stratification parameters of 
the external field-free plasma in the convection zone, with an extension of a sim- 
ple polytropic, sub-adiabatically stratified thin overshoot layer, as described by 
|Fan and Gong| (|2000p . The sub-adiabatic overshoot region extends from 
4.8xl0 10 cm (r = O.69i? ) to 5.026xI0 10 cm from Sun center, with the super- 
adiabatic convection zone extending from 5.026 xlO 10 cm to 6.75 xlO 10 cm (r = 
O.97i?0), which is still « 21 Mm below the solar surface. Profiles of T e , p e , p e , 
and the adiabaticity 5, we use can be found in Article 1. 

What sets the thin flux tube simulations in Article 1 (and continued here) 
apart from previous simulations is the inclusion into the drag force term of a 
time-dependent convective velocity field v e (r, t) relative to the rotating frame 
of reference. This three-dimensional global convection simulation is computed 
separately from the thin flux tube simulation using the ASH (anelastic spherical 
harmonic) code, as described by Micsch ~e£ aZ.I(|2006p (hereafter MBT06). In the 
anelastic approximation, the velocity of convection flows is taken to be much 
slower than the speed of sound in the fluid, and convective flows and thermal 
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variations are treated as a linear perturbation to a background state taken from 
a one-dimensional solar structure model. The computed v e (r,t) captures giant- 
cell convection and the associated mean flows such as meridional circulation and 
differential rotation, in a rotating convective envelope spanning r = O.69i?0 to 
r = O.97i?0 (4.8 xfO 10 cm to 6.75 xlO 10 cm from Sun center). ASH is a pseudo- 
spectral code, and the particular simulation used here is resolved by a grid of 
1 29 points in r, 256 points in 9, and 512 points in <f>. To find more specific details 
regarding this ASH simulation, see Article 1. 

Boundary conditions and simulation parameters for the ASH simulation are 
similar to Case AB3 in MBT06. Particularly, the lower thermal boundary condi- 
tion is the same as in Case AB3, with a latitudinal entropy gradient imposed to 
implicitly capture thermal coupling to the tachocline. This feature helps promote 
a conical rotation profile (see MBT06). The radial entropy gradient imposed at 
the outer boundary is steeper in this simulation, and therefore more in line with 
solar structure models (e.g. Christensen-Dalsgaard et aL[H996p . as compared to 
Case AB3. In this simulation the values of the thermal diffusivity K and the 
turbulent viscosity v at the outer boundary (r = 0.97i?©) are 4 x 10 13 cm 2 
s _1 and 2 x 10 13 respectively, and each decreases with depth in proportion to 
the inverse square root of the background density p~ x / 2 . Somewhat larger than 
the density contrast of Case AB3 due to a slightly deeper domain, the density 
contrast across this ASH simulation domain is 69, corresponding to 4.2 density 
scale heights. This yields a mid-convection zone Raylcigh number R a of 5 x 10 6 
and Reynolds number R e of order 50. The value of R e in this simulation is 
about a factor of two smaller than that used in the flux tube simulations of 
IJouve and Brunl (|2009[) . but the value of R a is larger. 

While this ASH simulation is more laminar than some others (|Miesch et alX 
2008; Jo uve and Brunl I2009|) . it still possesses all of the relevant features nec- 
essary to explore the fundamental interactions between thin flux tubes and the 
mean flows associated with global convection. These features include: asymmet- 
ric, rotationally aligned cells at low latitudes (density-stratified banana cells), 
rapidly-evolving downflows in the upper convection zone at high latitudes dom- 
inated by helical plumes, and a strong, solar-like differential rotation. Even the 
highest resolution simulations exhibit similar basic features, therefore we would 
not expect the essential results to change significantly with more turbulent 
convection. 

A typical giant-cell convection pattern at a depth of 25 Mm below the solar 
surface is shown in Figure [Tj Broad upflow cells are surrounded by narrow 
downflow lanes, which can reach maximum downflow speeds of nearly 600 m s _1 
at a mid-convection zone depth of about 86 Mm below the surface. Throughout 
the bulk of the convection zone, the combined influence of density stratification 
and the Coriolis force induces anti-cyclonic vorticity in expanding upflows and 
cyclonic vorticity in contracting downflows. These effects yield a mean kinetic 
hclicity density that is negative in the northern hemisphere and positive in 
the southern hemisphere (see Article 1 for an image of the associated kinetic 
hclicity). Such a helicity pattern is typical for rotating, compressible convection 
(e.g. IMiesch and Toomrel l2"009")l . 
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V r (m/s) 

Figure 1. A snapshot of the convective radial velocity at a depth of 25 Mm below the solar 
surface. 

Columnar, elongated downflow lanes align preferentially with the rotation axis 
at low latitudes, reflecting the presence of so-called "banana cells" (see Figurc[T|). 
Such structures propagate in a prograde direction relative to the polar regions, 
due in part to differential rotation and an intrinsic phase drift similiar to traveling 
Rossby waves (Micsch and Toomrel 12005]) . The total angular velocity fi/27r (with 
respect to the inertial frame) is solar-like, and decreases monotonically from 
«470 nHz at the equator to :=s330 nHz at the poles, and exhibits nearly conical 
contours at mid-latitudes (see Article 1 for an image of the associated differential 
rotation), as observed in the solar convection zone ( |Thompson et q/.]l2003p . 

As is in Article 1, our simulations start with isolated toroidal magnetic flux 
rings in mechanical equilibrium (neutral buoyancy) , located at a radial distance 
to the center of the Sun of r = ro = 5.05 x I 10 cm, slightly above the base of the 
convection zone. To ensure initial neutral buoyancy, the internal temperature of 
the flux tube is reduced compared to the external temperature. The toroidal 
ring is perturbed with small undular motions which consist of a superposition 
of Fourier modes with azimuthal order ranging from m — through m = 8 
with random phase relations. In this work and in Article I, the external time 
dependent three-dimensional convective flow described impacts the flux tube 
through its drag force term. We consider a range of initial field strengths of 15 
kG, 30 kG, 40 kG, 50 kG, 60 kG, and 100 kG, and initial latitudes ranging from 
1° to 40° for the toroidal flux ring, with a constant magnetic flux of 10 20 Mx to 
10 22 Mx. This range of magnetic flux values is typical of ephemeral regions and 
pores to the strongest sunspots. We also greatly improve statistics by performing 
seven simulations sampling different time ranges of the ASH convective flow for 
each magnetic field strength, initial latitude, and flux, for a total number of flux 
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tube simulations of 5940. The thin flux tube simulations sample a span of «2.7 
years of consecutive ASH convective flows. 

3. Results 

3.1. Rise Time and Latitude of Emergence 

The drag force term (the last term in Equation ((T|)) depends upon the magnetic 
field B and the flux $ through the ratio of (^/i?) 1 / 2 , which is proportional 
to the cross-sectional radius of the thin flux tube. As this ratio appears in the 
denominator of the drag force term, thinner tubes will experience more drag. 
For any given magnetic field, the diameter of the flux tube can be reduced 
by decreasing the magnetic flux. In the absence of convection, the drag force 
acting on the rising tube reduces the velocity of the flux tube in all directions, 
increasing its rise time and reducing its latitude of emergence. Figure [5} (a) shows 
a trend of increasing rise time of the flux tube simulations for both a decreasing 
magnetic field strength and a decreasing magnetic flux in the absence of con- 
vection, which has been found in previous thin flux tube simulations without 
the influence of convection {e.g. IMoreno-Insertisi 119831 IChoudhuri and Gilmanl 
IM71 ID'Silva and Choudhuril [TW31 IFan et al.\ 11555)1 . Taking these effects into 
consideration, a tube with a reduced flux and a reduced magnetic field should 
take the longest to rise, as is shown in Figure H]-(a), with 15 kG, 10 20 Mx tubes 
taking the longest time to rise. 

However, with the addition of convection, flux tubes with smaller magnetic 
flux exhibit a shorter rise time, except for those with an initial magnetic field 
strength of 100 kG, as shown in Figure [D(b). In this case, the drag force term 
affects how strongly the flux tube is coupled with the convective velocity field. 
Rising flux tubes with lower magnetic flux are advected strongly by convection. 
This aids the flux tube in emerging at the surface faster than it could at a 
larger magnetic flux, provided the drag force due to convection is significant 
compared to the buoyancy force. On the other hand, the magnetic buoyancy 
and magnetic tension of flux tubes with an initial magnetic field strength of 100 
kG dominate the drag force due to convection. As a result, flux tubes of 10 20 
Mx still take a longer time to emerge than those with a flux of 10 22 Mx, as 
is the trend without convection. Similar to Article 1, flux tubes with mid- field 
strengths of 40 kG - 50 kG and magnetic flux of 10 21 Mx and 10 22 Mx take the 
longest time to emerge. For these magnetic flux and magnetic field strengths, 
the average convective downflows and magnetic buoyancy of the flux tube are of 
similar magnitudes (see Article 1). A tug — of — war exists between these two 
effects until one eventually dominates. Considering the «1 1-year duration of the 
solar cycle, for flux tubes originating near the convection zone base, a maximum 
ruse time of about eight months for flux tubes subject to convective effects (see 
Fig. [2]-(b)) is much more realistic than a maximum rise time of about five years 
without convection (see Fig. H]-(a)). 

Without the external flow field, the increased drag force on flux tubes at lower 
magnetic flux aids in reducing the poleward deflection of the flux tube by the 
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Figure 2. (a) Average rise times for flux tube simulations without the influence of convection. 
Red crosses are for flux tubes with magnetic flux of 10 22 Mx, blue diamonds for 10 21 Mx, and 
black triangles for 10 20 Mx. Each symbol represents an average of nine flux tube simulations. 
This figure shows a decrease in emergence time with an increase in magnetic field strength and 
flux, (b) Average rise times for flux tube simulations with the influence of convection. Each 
symbol represents an average of Ri330 flux tube simulations. With the addition of convection, 
rise times of the flux tubes are shorter for tubes with smaller magnetic flux, except at 100 
kG where the trend is reversed. Flux tubes of mid-field strength now take the longest time to 
emerge. Bars represent the standard deviation of the mean. 



Coriolis force such that the tube rises more radially. As a result, the flux tube 
emerges at lower latitudes for lower magnetic flux, as has been found in previous 
thin flux tube calculations in the absence of convection 
(e.g. IChoudhuri and Gilmanl [19871 ID'Silva and Choudhuril IT9931 
IFan et ai[ I1993P . This is supported by Figure El which shows the latitudinal 
deflection (emergence latitude minus initial latitude) of the flux tube as a func- 
tion of initial latitude in the absence of convection for magnetic flux of 10 20 Mx 
(black triangles), 10 21 Mx (blue diamonds), and 10 22 Mx (red crosses). Also, as 
found in previous studies (IChoudhuri and Gilmanl 119871 |Caligari et al.\ I1995[) , a 
low-latitude zone with the absence of flux emergence due to poleward deflection 
is found for magnetic field strengths of 15 kG to 40 kG. As such, these flux tubes 
would not be able to produce active regions near the equator. 

With convection, the previous problem of poleward slippage for flux tubes 
of weak initial field strength is rectified. Convection produces a scatter in the 
emerging latitude compared to the cases without convection, which is both pole- 
ward and equatorward. Figure 0] shows the latitudinal deflection of the flux tube 
for fluxes of 10 20 Mx, 10 21 Mx, and 10 22 Mx resulting from simulations with 
convection (plus signs), compared to the case without convection (diamonds). 
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Figure 3. Latitudinal deflection (emergence latitude minus initial latitude) of the flux tube 
apex as a function of initial latitude for initial flux tube magnetic field strengths of 100 kG, 
60 kG, 50 kG, 40 kG, 30 kG, and 15 kG in the absence of convection for magnetic flux of 10 22 
Mx (red crosses), 10 21 Mx (blue diamonds), and 10 20 Mx (black triangles). An increased drag 
force for tubes with a smaller flux reduces the poleward deflection of the tube. 



A negative latitudinal deflection indicates that the flux tube emerges close to 
the equator than where it originally started. Such a phenomenon only occurs 
due to external convective flows. It is even possible for portions of the flux tube 
to cross into the opposite hemisphere. We note a systematic trend of reduced 
poleward deflection of the flux tube for all magnetic fluxes considered here. 
This is in agreement with Article 1, where we only studied 10 22 Mx cases. It is 
also evident that there is a greater amount of scatter in latitudinal deflection 
at lower magnetic flux, which occurs because the tube is advected more by 
convective flows than at a larger magnetic flux. Also, as the magnetic field 
strength increases for a particular flux, the latitudinal deflection decreases, as 
is the case for tubes both with and without convection. This is a result of the 
buoyancy force overpowering the Coriolis force at large magnetic fields, forcing 
the tube to rise more radially. 
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Figure 4. Latitudinal deflection (emergence latitude minus initial latitude) of the flux tube 
apex as a function of initial latitude for initial flux tube magnetic field strengths of 100 kG, 
60 kG, 50 kG, 40 kG, 30 kG, and 15 kG with a magnetic flux of (a) 10 20 Mx in column 1, (b) 
10 21 Mx in column 2, (c) and 10 22 Mx in column 3. Red diamond are for flux tubes without 
convection, and plus signs arc for flux tubes with convective effects. Both axes are in units of 
degrees. Addition of a convective velocity field results in flux tubes that are able to emerge 
near the equator at lower magnetic fields. 
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3.2. Joy's Law 

The phenomenon associated with solar active regions known as Joy's Law de- 
scribes the tilting behavior of emerging flux regions toward the equator (e.g. 
Hal e et aL[ ITT) 1 Pi. A line drawn between the center of the two bipolar regions 
will be tilted with respect to the east- west direction such that the leading polarity 
(in the direction of solar rotation) is closer to the equator than the following. On 
average, the tilt of this line with respect to the East -West direction, quantified 
by the tilt angle, increases with an active region's latitude of emergence. In this 
article, the tilt angle is computed as the angle between the tangent vector at 
the apex of the emerging loop (once it has reached the top of the simulation 
domain), and the local East -West direction. We define a positive sign of tilt 
as a clockwise (counter-clockwise) rotation of the tangent vector away from the 
East -West direction in the northern (southern) hemisphere, consistent with the 
direction of the observed mean tilt of active regions. If the magnitude of the 
tilt angle exceeds 90°, then the active region is of an anti-Hale arrangement, 
possessing the wrong leading polarity in the direction of solar rotation for that 
particular hemisphere. This is the same convention as used in Article 1. 

Tilt angles of the flux tube as a function of emergence latitude for magnetic 
flux of 10 20 Mx, 10 21 Mx, and 10 22 Mx are shown in Figure[SJ for initial magnetic 
field strengths of 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG both with 
(plus signs) and without convection (diamonds) . Each of the panels in this figure 
includes about 3.5 times more flux tube simulations per magnetic field strength 
compared to Article 1, sampling different time ranges of the ASH convective 
velocity flow, therefore greatly improving statistics. As a result, each panel 
contains «330 points. Here we also expand the investigation to include cases 
with flux of 10 20 Mx and 10 21 Mx. As in Article 1, we assume that the tilt angle 
increases monotonically with increasing emergence latitude, and so perform a 
linear least-squares fit of Joy's Law to the data using the equation: a — toaA, 
where a, m^, and A represent the tilt angle, slope, and latitude respectively. The 
fit is forced to go through zero because no tilt is expected for equatorial sunspot 
groups. Figure [5] shows in gray the best- fit line for the data without convection, 
and black for the flux tubes with convection. The slopes of these best-fit lines 
along with their uncertainties are reported in Table [TJ for flux tubes without 
convective effects, and Table [2] for flux tubes with convective effects. 

Without the influence of convection, untwisted flux tubes tilt in the appro- 
priate direction for their respective hemispheres due to the Coriolis force acting 
on the limbs of the emerging loop (e.g. ID'Silva and Choudhuril ITTT9T3"]) . However, 
for weaker initial field strengths of 15 kG and 30 kG, some emerging loops show 
negative tilt angles, opposite to the sign of the active region mean tilts. This 
occurs because plasma flow along the flux tube near the apex changes from di- 
verging to converging as it enters the upper convection zone (e.g. |Caligari et al\ 
1995] IFan and Fisherl I1996p . The Coriolis force acting on the converging flow 
drives a tilt of the wrong sign, resulting in a negative best fit line slope for tubes 
without convection at magnetic field strengths of 15 kG and 30 kG. This effect 
is especially severe here for flux tubes of 15 kG and 10 20 Mx. The tilt angle of 
these tubes is so large and of the wrong sign for initial latitudes of 5°, 10°, and 
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Figure 5. Tilt angles as a function of emergence latitude for initial magnetic field strengths 
of 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG for cases with (plus signs) and without 
(diamond points) the influence of convection for magnetic flux of (a) 10 20 Mx in column 1, (b) 
10 21 Mx in column 2, and (c) 10 22 Mx in column 3, with data sets sampling seven different 
convective velocity flow fields. The gray line is the linear best fit for flux tubes in the absence 
of convection, and the black line is the best-fit line for the cases subjected to convective flows. 
Slopes for all of these best- fit lines are shown in Table 1. A color bar indicates the original 
starting latitude of the flux tube. All axes are in units of degrees. Convection introduces a 
scatter of the tilt angle about the best-fit line, and also aids in increasing the slope of the 
best-fit line, especially at lower magnetic field strengths. 
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Table 1. Slopes uia of the linear best-fit lines to the tilt angle as 
a function of emergence latitude for flux tubes without convective 
effects. At low magnetic field strengths, the slope of the best-fit line 
is negative, indicating a departure from the Joy's Law trend. 



B 10 20 Mx 10 21 Mx 10 22 Mx 

without convection 



100 kG 


0.10 


± 


0.01 


0.16 


± 


0.01 


0.25 


± 


0.01 


60 kG 


0.18 


± 


0.01 


0.27 


± 


0.01 


0.30 


± 


0.02 


50 kG 


0.18 


± 


0.01 


0.25 


± 


0.01 


0.22 


± 


0.02 


40 kG 


0.16 


± 


0.01 


0.16 


± 


0.01 


0.12 


± 


0.01 


30 kG 


-0.12 


± 


0.01 


-0.08 


± 


0.03 


-0.07 


± 


0.05 


15 kG 


-0.50 


± 


0.38 


-0.24 


± 


0.22 


-0.13 


± 


0.13 



Table 2. Slopes mA of the linear best-fit lines to the tilt angle 
as a function of emergence latitude for flux tubes with convec- 
tive effects. Convection aids in increasing the slope of the best-fit 
line, and the slope also peaks at mid-field strengths of 40 kG 
and 50 kG. In comparison, using white light sunspot group tilt 
angle data, Dasi-Espuig et al. (2010) find ra\ = 0.26 ± 0.05 for 
Mount Wilson data and m . 4 = 0.28 ± 0.06 for Kodikanal data. 



B 10 20 Mx 10 21 Mx 10 22 Mx 

with convection 



100 kG 


0.36 


± 


0.05 


0.28 


± 


0.02 


0.25 


± 


0.02 


60 kG 


0.44 


± 


0.08 


0.37 


± 


0.04 


0.35 


± 


0.02 


50 kG 


0.50 


± 


0.08 


0.38 


± 


0.04 


0.42 


± 


0.04 


40 kG 


0.36 


± 


0.10 


0.45 


± 


0.08 


0.44 


± 


0.06 


30 kG 


0.44 


± 


0.11 


0.31 


± 


0.11 


0.37 


± 


0.09 


15 kG 


0.29 


± 


0.12 


0.18 


± 


0.14 


0.27 


± 


0.12 



15°, that they do not appear on the tilt angle plot (Figure 0(a), lower left). 
As can be seen from Table [1] and Table [2j we find that at all magnetic fields, 
the slopes of the best-fit lines for the flux tubes with convection are increased 
from the cases without convection. This occurs in part because the upflows in 
the convective velocity field have an associated kinetic hclicity that helps to 
drive the tilt of the flux tube apex in the appropriate Joy's Law direction for 
its respective hemisphere (i.e. toward the equator). Such a mean kinetic hclicity 
corresponds to a vertical vorticity in upflows that is clockwise in the northern 
hemisphere, and counter-clockwise in the southern hemisphere. 

With convection and for all magnetic fluxes, all magnetic field strengths of 15 
kG - 100 kG show positive best-fit line slopes, agreeing with Joy's Law. This is 
an improvement upon Article 1, for which there were fewer simulations only for 
a magnetic flux of 10 22 Mx and the best-fit line slopes for 15 kG and 30 kG of 
0.12 ± 0.18 and 0.15 ± 0.21 respectively, had too large uncertainties to report 
a definitive Joy's Law trend. It is interesting to note that with convection, the 
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mid-magnetic field strength flux tubes (40 kG - 50 kG) have the largest best-fit 
slopes. This is probably due to the fact that at these mid-field strengths, the 
joint effects of magnetic buoyancy and the convective flows are such that the 
tubes have a longer rise time compared to the other field strengths. Therefore, 
the systematic effects from the Coriolis force and kinetic helicity in convective 
upflows have a longer time to act on the flux tubes, and the flux tube will likely 
encounter multiple convective cells during its evolution. Since the mean kinetic 
helicity is obtained by averaging the kinetic helicity of many cells over time, a 
flux tube which takes a longer time to emerge will be influenced more by the 
mean kinetic helicity rather than stochastic fluctuations. 

For solar observational data, the observed tilt angle of the active region is 
often plotted as a function of its emergence latitude for active regions of all sizes 
and magnetic flux. Figure [5] shows the tilt angles from the flux tube simulations 
together for selected values of magnetic field and flux. A linear least-squares fit 
is performed on this data for various flux and magnetic field combinations as 
shown in Figure [6j using the fit of a — m^A. From white light sunspot group 
tilt angle data spanning solar cycles 15 - 21, |Dasi-Espuig et aL\ (|2010|) find a 
linear best-fit line slope of 0.26 ± 0.05 for Mount Wilson data, and 0.28 ± 
0.06 for Kodikanal data. For simulated flux tubes of all magnetic field strengths 
considered and fluxes of 10 21 Mx and 10 22 Mx (Figure [6]- (b)), as these are most 
likely the magnetic flux strengths required to produce sunspot groups that can 
be identified in white light, we find a slope of = 0.34 ± 0.02, which docs 
fall within the range of tilt angles that |Dasi-Espuig et ai\ (|2010p find given the 
uncertainties on the slope. 

A Joy's Law fit can also be performed using a = mssin(A), which is a good 
choice assuming the origin of the tilt angle is related to the Coriolis force, as 
this force varies with latitude as sin(A). Here, raj is the best-fit line slope. 
IStenflo and K osovichevl (|2012[) perform such a fit, which gives a slope of mj = 
32.1° ± 0.7°, using 15 years of MDI full-disk magnetograms. Wang and Sheele"y| 
(1983|) find a similar trend using National Solar Observatory Kitt Peak mag- 
netograms of bipolar active regions for Solar Cycle 21. Performing the same fit 
on the tilt angles of our simulated flux tubes for magnetic fluxs of 10 20 Mx, 
10 21 Mx, and 10 22 Mx (Figure [6} (c)), which encompasses the range of flux these 
magnetograms consider, we find raj = 22° ± 1°. This value is much lower than 
what |Wang and Sheeley] (|1989[) and IStenflo and Kosovichevl (|2012p determine. 
However, using Mount Wilson white light sunspot group data collected from 
1917 - 1985. IFisherlTlID ([1995)) find m B = 15.69° ± 0.66°. For our simulated 
flux tubes of 10 21 Mx and 10 22 Mx (Figure [i-(b)), we find m B = 21° ± 1°, which 
is slightly larger than the value derived by Fis her et al.\ (|1995p . Recalling that 
mid- field strengths of 40 kG - 50 kG have the largest linear best-fit slopes m^, 
we perform the fit a = mssin(A) only for these field strengths and all magnetic 
fluxes we consider (Figure EJ-(d)). The result is a significantly increased value of 
m B = 26° ± 2°, although it still falls short of the 32.1° ± 0.7° value found by 
IStenflo and Kosovichevl (|2012p . For all values of magnetic flux and magnetic field 
that we consider for Figures [6} (a) through EJ-(d), the best- fit lines we obtain for 
our simulated flux tubes arc bounded by the tilt angle best-fit lines as derived 
from observations. We note that the best-fit lines for the Joy's Law tilt angle 
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dependence derived from white light sunspot group data (jFisher et all [1995; 



Dasi-Espuig et al. I2010p are much lower than those derived from magnetograms 



QWang and Sheeley[ UM§[ IStenflo and Kosovichevl I20T2|1 . This night occur be- 



cause weaker active regions can be identified in magnetograms, which may not 
appear in white light sunspot group images. However, it may also be the result of 
selection effects employed by various authors, as well as the fact that anti-Hale 
regions might not be identified from white light data. 

According to IFan et ai\ fjl994[) , in the absence of convection, the tilt angle 
of the flux tube should increase with increasing flux as a result of the Coriolis 
force acting on the buoyantly rising flux tube. With convection, given the un- 
certainties that we obtain, we do not find a statistically significant dependence 
of the Joy's Law slope on magnetic flux (see Tabled]), which agrees with the 
results of IStenflo and Kosovichevl (|2012p . We have shown that the addition of 
convection increases the Joy's Law slope, regardless of the amount of time that 
the flux tubes stays in the bulk of the convection zone. This indicates that 
convective effects can have a significant contribution to the tilt angle of the flux 
tube throughout its evolution. Therefore, the results of IStenflo and Kosovichevl 
(2012]) do not rule out the paradigm that flux tubes obtain at least a portion of 
their tilt angle during their buoyant rise. 

We also note that 6.9% of our simulated flux tubes (see Figure 0(c)) emerge 
with polarites that violate Hale's Law, in comparison to the «4% as found 
via observations of medium to large-sized active regions ( |Wang and Sheeley] 
119891 IStenflo and Kosovichevl 120120 . In our simulations, these violations arise 
as a result of flux tubes emerging in the opposite hemisphere from which they 
originated, or as a result of the flux tube becoming so distorted by convection 
that the legs of the emerging loop can become reversed. Both of these situations 
usually occur for weak magnetic field strength cases. The majority of emerging 
anti-Hale flux tubes with moderate initial magnetic fields between 40 kG-50 
kG happen as a result of hemisphere crossing due to convective flows. Therefore, 
in light of our findings, observations of anti-Hale spots of various magnetic flux 
does not rule out the rising flux tube paradigm. 

3.3. Tilt Angle Scatter 

We attempt to further constrain the magnetic field strength at which the solar 
dynamo might be operating by observing the scatter of tilt angles our simulated 
flux tubes produce about the Joy's Law trend. In Figure [5J there is a clear 
scatter of the tilt angle about the best-fit line for the data with convection, 
which becomes greater at lower magnetic flux. Again, this trend is the result 
of stronger coupling between the flux tube and convective flows at a reduced 
magnetic flux. To quantify the scatter of the tilt angles around the best fit line, 
we calculate the standard deviation of the tilt about its fitted value: 



(6) 
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Figure 6. Tilt angles for various flux and magnetic fields plotted together. The dark blue 
line is the linear best-fit line for tilt angle as a function of latitude as found from white light 
sunspot group observations (Dasi-Espuig et al., 2010), and the green line is the best-fit line 
for tilt angle as a function of the sine of latitude as found from magnctograms (Stenflo and 
Kosovichev, 2012). In all graphs, we also plot the best- fit line of the simulations for tilt angle 
as a function of latitude (red, slope A), and for the tilt angle as a function of the sine of 
latitude (light blue, slope B) for (a) tilt angles of all magnetic field strengths and flux of 10 20 
Mx only, (b) tilt angles for all magnetic field strengths and fluxes of 10 21 Mx and 10 22 Mx, 

(c) tilt angles for all magnetic field strengths and fluxes of 10 20 Mx, 10 21 Mx, and 10 22 Mx, 

(d) tilt angles for all fluxes for magnetic field strengths of 40 kG-50 kG. The best-fit lines for 
the thin flux tube simulations are bounded by the observational best-fit lines. We find that 
6.9% of our simulated flux tubes are of an anti-Hale configuration, as compared to the Ri4% 
as found by Stenflo and Kosovichev (2012). 
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Table 3. Standard deviation of the tilt angle about its fitted value. The standard deviation 
tends to increase with decreasing flux. 





15 kG 


30 kG 


40 kG 


50 kG 


60 kG 


100 kG 


10 20 Mx 


63 ± 4 


58 ± 3 


54 ± 3 


41 ± 2 


40 ± 2 


22 ± 1 


10 21 Mx 


73 ± 4 


55 ± 3 


42 ± 2 


30 ± 2 


19 ± 1 


10.4 ± 0.6 


10 22 Mx 


71 ± 4 


52 ± 3 


35 ± 2 


23 ± 1 


12.5 ± 0.7 


8.2 ± 0.4 


10 21 , 10 22 Mx 


74 ± 2 


54 ± 2 


39 ± 2 


26 ± 1 


16.3 ± 0.6 


9.4 ± 0.4 


10 20 , 10 21 , 10 22 Mx 


69 ± 2 


55 ± 2 


44 ± 1 


32 ± 1 


26.4 ± 0.8 


14.9 ± 0.5 



where a>i is the i tilt angle, a/u is the i tilt angle as a result of the fit for 
which the slopes ttla are reported in Table [H and N is the number of points 
considered. 

Using Mount Wilson white light sunspot group data (1917-1985), 
IFisher et al\ (TIMS)) find that a < 40° for all emergence latitudes. We evaluate 
a for each field strength and flux value, and the results are given in Tabic [21 
Comparing to the results of IFisher et a l. (1995), this would exclude magnetic 
field strengths of 15 kG and 30 kG for the combination of fluxes of 10 21 Mx and 
10 22 Mx together (row 4 in Table [3]) as the progenitors of solar active regions, 
which has a scatter that exceeds the observed value. We only consider these flux 
values because IFisher et al.\ (|1995p uses white light sunspot group data. If we 
include all magnetic fields and these two fluxes together (not shown in Table 
we obtain a a value of 41.4 ± 0.7, which is fairly close to what lFrsher et aL\ 
(1995) obtain. However, if we exclude magnetic fields of 15 kG and 30 kG for 
fluxes of 10 21 Mx and 10 22 Mx (not shown in Tabic [3]), wc obtain a a value of 
25.3 ± 0.5, well within the range that IFisher et aZ.1 (|1995j> suggests. 

Active regions of flux values of 10 20 Mx would appear on magnctograms, but 
would probably not be associated with white light sunspot groups. Since we 
have simulations for 10 20 Mx, we calculate a for these values as well in Table 
[3] For this magnetic flux, a is very large and exceeds the values suggested by 
IFisher et all ©95) for all magnetic field strengths <40 kG. However, this result 
does not exclude these magnetic field values as progenitors of active regions 
for a flux of 10 20 Mx, since these regions most likely would not appear in 
the IFisher et al.\ (|1995[) study. We do note that while the Joy's Law trend is 
not statistically dependent on flux (see Tabled, the scatter of the tilt angles 
about their best-fit line tends to increase with decreasing flux (see Table [3]). 
|Wang and Sheeley| (TT9"g9")) and lStenflo and Kosovichevl l|20T2)l also find that the 
spread in the distribution of the tilt angles increases as one goes to smaller 
regions, or smaller magnetic flux. 

3.4. Preferred Tilt Angle 

IHowardl (|1996[) investigated the distribution of sunspot group tilt angles from a 
set of Mount Wilson white light photographs (1917-1985) as well as sunspot 
groups derived from plage in Mount Wilson daily magnctograms (1967-1995). 
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Figure 7. Distribution of tilt angles in 2.5° bins for all initial magnetic field strengths and 
all magnetic fluxes considered in this study. The peak of the distribution function is between 
10° - 12.5° (black line), with a non-linear least-squares Gaussian fit (red line) of center 10.6°, 
as shown by the dash-dotted line. 

They bin the tilt angles in 2.5° increments, and plot the number of ocurrances 
of the tilt angles in these bins. In this way, they find that sunspot groups tend 
to form most often with tilt angles between 2.5° -5° for white light data, and 
between 7.5° -10° for magnetogram data. The tilt angle distribution that we 
obtain for flux tubes of all magnetic fields and magnetic flux studied here using 
tilt angle bins of 2.5° is shown in Figure [3 We find that for all magnetic fluxes 
considered, the tilt angle distribution peaks between bins of 10° - 12.5° (black 
line), where a non- linear least-squares Gaussian fit (red line) gives a center, or 
preferred tilt, of 10.6°. This is in good agreement with the magnetogram data 
of IHowardl (|1996p and IStenflo and Kosovichevl (|2012p . who find a preferred, or 
most common, tilt of «10° for all bipolar regions emerging between a 15° -20° 
latitude range including the largest active regions down to regions of « 10 20 Mx. 

IHowardl (T1996) also finds an average tilt angle of 4.3° for sunspot groups and 
6.3° for the magnetograms, whereas we find an average of 9.5° ± 0.6° for all 
fluxes considered together, where the uncertainty is the standard deviation of 
the mean. The average tilt angle value is significantly larger than what IHowardl 
(199(5]) finds. However, this could be due to the nature of our chosen convective 
velocity simulations, the thin flux tube idealization, or a result of the selection 
criteria that IHowardl (|1996[) uses to identify sunspot groups or plages. They 
count an active region on every day it is observed, so an active region may be 
represented a number of times for a single solar disk passage. As such, their data 
also includes the evolution of the region's tilt angle as it ages. 

3.5. Magnetic Fields at the Flux Tube Apex 

Using thin flux tube simulations, IFan et al\ (|1993p show that the proceeding 
leg of an emerging flux loop has a stronger magnetic field than the following 
leg as a result of the Coriolis force, which induces a differential stretching of 
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the rising loop. This provides an explanation for the observed coherent, less 
fragmented morphology for the leading polarity flux of an active region (e.g. 
|Bray and Lough head I1979p . Subsequent simulations using the more physical 



mechanical equilibrium state (e.g. Caligari et al. 119951 IFan and Fisheri [T996). 



as opposed to temperature equilibrium (jFan et aU 11993ft . found that the leading 
leg of the emerging loop has a stronger magnetic field than the following only 
for flux tubes with an initial field strength below 60 kG. As is in Article 1, we 
investigate the magnetic field asymmetry by calculating dB/ds, the derivative 
of the magnetic field along the arc length s in the direction of solar rotation, 
at the apex of the emerging loop. If dB/ds is greater (less) than zero, then 
the leading (following) leg has a stronger magnetic field. We find that for all 
fluxes, tubes with an initial magnetic field < 50 kG tend to emerge with the 
appropriate magnetic field asymmetry of a stronger field strength in the leading 
side, consistent with the tendency of a more coherent leading polarity in an 
active region. 

Recent three-dimensional radiation MHD simulations of magnetic flux emer- 
gence and sunspot formation in the top layer of the convection zone and the 
photosphere by Rcmpel (private communication, 2012) have suggested that the 
morphological asymmetry of sunspots may be caused by a flow of plasma along 
the emerging tube out of the leading portion of the active region into the follow- 
ing leg, and is less dependent on the asymmetry in field strength of the emerging 
loop. If this is indeed the case, then we can not rule out flux tubes with stronger 
initial magnetic field strengths of > 50 kG at the base of the solar convection 
zone as the progenitors of solar active regions, simply due to their magnetic field 
asymmetries. 

The nature of the magnetic field below the visible surface of the Sun is not 
well known. Simulations of flux emergence in the upper layers of the convection 
zone have been computed in domains that span from the photosphere to depths 



of 2 Mm to 20 Mm below the surface (e.g. Cheung et al. 120071 Rempel et al, 
2009; St ein et aDl2011|) . The simulations assume a range of initial magnetic field 
strengths of the flux tube at the bottom of their simulation domain anywhere 
from «2 kG to «30 kG. We show the magnetic field strength at the apex of our 
thin flux tube simulations, once the tube has reached the top of the simulation 
domain at a depth of 25 Mm below the solar surface, for multiple initial magnetic 
field strengths and flux in Figure [8] Our results do support average magnetic 
field strength values of « 500 G to w 15 kG for flux tubes reaching a height of 
«21 Mm below the solar surface, assuming initial field strengths of 15 kG to 100 
kG at the base of the convection zone for magnetic fluxes of 10 20 Mx to 10 22 
Mx. The magnetic field strength of the flux tube at the apex is largest for larger 
initial fields, and the scatter of the magnetic field at the apex for a particular 
magnetic field strength increases with decreasing flux. 

3.6. Emerging Loop Rotation Rate 

Observations show that sunspots tend to rotate faster than the solar surface 
plasma ( Howard and Harvey, 1970| |Golub and Vaiana, 1978 ). The azimuthal 
velocity of leader, follower, and all sunspots (in the reference frame with a solid 
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Figure 8. Magnetic field strength at the apex of the flux loop for (a) 10 20 Mx in column 1, 
(b) 10 21 Mx in column 2, (c) and 10 22 Mx in column 3. Our results support average magnetic 
field strengths of f»500 G to Kil5 kG for flux tubes apices which have reached the top of the 
simulation domain, which is still 25 Mm below the solar surface. 
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Figure 9. Rotation rate of leader, follower, and all sunspots, plotted with the observed az- 
imuthal rotation rate of the solar plasma at the surface (blue dash-dotted line) and r = O.95R0 
(red dashed line). All values are plotted with reference to the solid body rotation rate of 
f2o=2-7 X 10 — 6 rad s -1 , so that the zero line is the solid body rotation of the Sun. This image 
suggests that sunspots rotate at nearly the same rate as the solar plasma at O.95-R0. 



body rotation rate of Slo = 2.7 x 10 -6 rad s _1 ) as derived from 
IGilman and Howardl (|1985p are plotted in Figure [HJ as well as the observed az- 
imuthal rotation rate of the solar plasma at the surface (blue line) as determined 
from surface spectroscopic Dopplcr-vclocity measurements 



(e.g. Thompson et al. I2003|) and the rotation rate at r = O.95i?0 (red line) as 
found via inversions of helioseismic observations ( |Howe et al., 2000| ) . This image 
suggests that sunspots tend to rotate at nearly the same rate as the solar plasma 
at O.95i?0, and therefore rotate faster than solar surface plasma. This might also 
indicate that sunspots are anchored at a depth of O.95i?0, near the surface shear 
layer interface. 

In order to quantify the rotation rate of emerging loops from our thin flux 
tube simulation, we calculate the apparent azimuthal speed of the center point 
between the leading and following intersections of the emerging loop with the 
constant r surface of 0.95i?Q during the last two time steps before the loop apex 
reaches the top of the simulation domain. This reflects the apparent azimuthal 
speed of an emerging active region, and it is not the actual azimuthal speed of the 
flux tube plasma at the tube apex. We do this only for flux tubes of 10 21 Mx and 
10 22 Mx, assuming tubes of flux order 10 20 Mx do not produce large sunspots. 
|Caligari et al. \ (|1995p computed the azimuthal phase speed of the summit of an 
emerging flux loop without convective effects throughout its rise. This phase 
motion is related to the wave character of the rising flux loop, which is similar 
in behavior to a transversal wave propagating along a string. They find that the 
azimuthal phase speed decreases with increasing height of the summit, and that 
it changes sign at about 50 Mm below the solar surface such that the angular 
velocity of the summit is smaller than that of the external plasma. However, 
they suggest that the resulting active region will still show a higher rotation 
velocity than the surrounding plasma because of the inclination difference of the 
leading and following legs with respect to the local vertical. This will occur when 
the following leg has a steeper slope than the leading leg, which is caused by 
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the conservation of angular momentum as the tube rises ([Moreno- Insertis et al.i 
IT9"M1 |Caligari et at] UMB |Caligari et at] HM^i . In Article 1, we find that, on 
average, all magnetic field strengths between 15 kG and 100 kG at a magnetic 
flux of 10 22 Mx show the appropriate inclination asymmetry of the emerging flux 
loop such that the leading and following polarities of a resulting active region 
will appear to drift apart from each other and rotate faster than the surrounding 
plasma. To make comparisons to actual sunspot rotation rates, we feel that the 
method presented in this section is a better alternative. 

Figures [TO]- (a) andfTDlfb) show the average of the apparent azimuthal speed in 
5° bins, for cases with a magnetic flux of 10 21 Mx and 10 22 Mx respectively, where 
bars on the points are the standard deviation of the mean. For comparison, these 
plots also show the average azimuthal rotation rate of the ASH simulation at r = 
O.95i?0 (red line), and the azimuthal flow speed one would expect at the surface 
(blue line) assuming the differential rotation in the top shear layer between the 
surface and 0.95-R© decreases by «10 nHz (what we will subsequently call the 



inferred surface rate) as found by hclioseismology (e.g. Thompson et al. I2003|) . 

In all cases, flux tubes that emerge at high latitudes above 40° rotate at a rate 
that is close to, or faster than, the ASH rate at r = O.95i?0 and hence faster than 
the inferred surface rate expected of our simulation. With initial magnetic field 
strengths of 15 kG and 30 kG, for both magnetic fluxes, the average rotation 
rates of the majority of flux tubes that emerge at < 35° are less than the inferred 
surface rate. However, at mid-magnctic-field strengths of 40 kG and 50 kG, the 
rotation rates of the flux tubes roughly follow the inferred surface rate. Only for 
initial magnetic field strengths of > 60 kG are the flux tubes capable of rotating 
at or faster than the inferred surface rate for all emergence latitudes, considering 
magnetic fluxes of 10 21 Mx and 10 22 Mx. 

We recognize that this method of quantifying the rotation rate of the emerging 
flux tubes in our simulation has some limitations. Due to the nature of the ASH 
convection simulation and the fact that the thin flux tube simulation breaks 
down in the upper portion of the convection zone, it is not possible to allow the 
thin flux tube to emerge all the way to the solar surface. Instead, we operate 
under the assumption that the rotation rate of the emerging loop at a constant 
r surface of O.O5i?0 below the solar surface is a good representation for how 
an active region will behave at the solar surface, as is reflected in Figure [H] for 
solar observations. The discrepancy between the observed and simulated rotation 
rates could in part be attributed to the particular ASH simulation that we use, 
which does not precisely reproduce the solar f2-profile. Our simulated flux tubes 
are effectively anchored at the convection zone base, whereas real flux tubes may 
decouple from the deep convection zone at some point during their evolution, 
and become anchored closer to the surface. This is why we choose to investigate 
the rotation rate closer to the surface. As our thin flux tube is a one-dimensional 
string of mass elements, we cannot address the issue that the flux tube could 
lose its coherency and become fragmented ( |Longcope et al.\ 11996ft , which may 
result in a stronger coupling between the tube and convective fluid motions, and 
could be a significant contributing factor to the rotation rate of active regions. 

In an attempt to gain a stronger coupling between the flux tube and convective 
fluid motions, we perform some simulations where the drag coefficient Cd in the 
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Figure 10. Average rotation rate of emerging flux loops subjected to the convective flow 
with magnetic flux of (a) 10 21 Mx on the right (rightmost six panels) and (b) 10 22 Mx on the 
left (leftmost six panels). Averages are taken in 5° bins, with bars representing the standard 
deviation of the mean. Red and blue lines correspond to the ASH rotation rate at r = 0.95-R© 
and inferred solar surface rotation rate respectively, assuming a difference of 10 nHz between 
the two. All values are plotted with reference to the solid body rotation rate of Qo = 2.7 X 10 — 6 
rad s , so that the zero line is the solid body rotation of the Sun. Strong flux tubes with 
initial field strengths of > 60 kG rotate faster than or nearly equal to the inferred surface rate. 



last term of Equation ([T]) is increased from unity to the constant values of 1.5 
and 2. Also, we perform simulations where we adjust the drag coefficient such 
that it exponentially varies from 1 at the base of the convection zone to 1.5 
or 2 at the top of the simulation domain. These tests were performed for flux 
tubes with initial field strengths of 40 kG and magnetic flux of 10 22 Mx. These 
efforts to alter the drag coefficient to produce stronger coupling to convection did 
not result in faster rotation rates compared to the surface rate for all latitudes 
less than 25°. This indicates that 40 kG flux tubes or less will probably not be 
able to reproduce sunspot rotation rates utilizing the thin flux tube model as it 
currently stands. 



4. Summary 

By embedding the thin flux tube model in a three-dimensional, turbulent, con- 
vective velocity field representing the solar convective envelope, we study how 
convection can influence the properties of emerging active region flux tubes. In 
comparing these properties to those obtained from solar active region observa- 
tions, we attempt to constrain the magnetic field strengths of dynamo-generated 
magnetic fields at the base of the solar convection zone. The thin flux tube ap- 
proximation, although idealized, allows us to investigate active region scale flux 
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tubes at weak to strong magnetic field strengths of 15 kG- 100 kG under perfect 
frozcn-in flux conditions. Wc find that subjecting the thin flux tube to turbulent 
convective flows does indeed alter flux tube dynamics, and that it can have a 
significant impact on the properties of the emerging flux loop in comparison to 
flux tube simulations performed in the absence of a convective velocity field. 
Also, the addition of convection aids the flux tube in more closely reproducing 
some observed properties of solar active regions. This article extends upon the 
previous work of Article 1 by including flux tube simulations with additional 
magnetic fluxes of 10 20 Mx, 10 21 Mx, as well as 10 22 Mx, and increasing the 
number of simulations performed per magnetic field strength by a factor of 3.5 
in order to improve statistics such that uncertainties are reduced. We also include 
more observational diagnostics, such as tilt angle scatter and sunspot rotation 
rate, to put further constrains on the field strength of the dynamo-generated 
magnetic field as the progenitors of solar active regions. 

Decreasing the magnetic flux of the tube results in an increase of the drag 
force acting on the rising flux loop. With convection, flux tube rise times decrease 
with decreasing flux for initial magnetic field strengths of 15 kG - 60 kG because 
the increased drag force causes the flux tube to be more closely coupled with 
convection, and so they are advected by turbulent flows more strongly than the 
10 22 Mx case. For all magnetic fluxes that we consider here, flux tubes are able 
to emerge near the equator with the aid of convective flows, which solves the 
previous problem of poleward slip for flux tubes of low magnetic field strengths 
without convective effects. 

With the increased number of thin flux tube simulations, and hence improved 
statistics in this study, we are able to confidently say that for all magnetic field 
strengths of 15 kG to 100 kG, and all magnetic fluxes studied here, produce 
emerging loops with tilt angles that follow the Joy's Law trend. This is an 
improvement upon Article 1 where the slope of the linear Joy's Law fit of the 
tilt angles of the emerging loops as a function of emerging latitude had too large 
an uncertainty for magnetic field strengths of 15 kG and 30 kG to report a 
definitive Joy's Law trend. Of particular note is the fact that helical convective 
upflows help to drive the tilt angle of the flux tube in the appropriate Joy's law 
direction for both hemispheres. Including all tilt angles together for all magnetic 
field strengths and magnetic flux of 10 21 Mx and 10 22 Mx, we calculate the linear 
best-fit line for the tilt angle as a function of emergence latitude, and find a slope 
of m A = 0.34 ± 0.02, which overlaps with the range of 0.26 ± 0.05 and 0.28 ± 
0.06 as suggested by white light sunspot group image from Mount Wilson and 
Kodikanal, respectively ( |Dasi-Espuig et oT|l2010p . Performing a fit for the tilt 
angle as a function of sine of the latitude for 10 20 Mx, 10 21 Mx, and 10 22 Mx, 
we find a best-fit line slope of ms — 22° ± 1°, which is greater than the value 
15.69° ± 0.66° obtained from white light sunspot group data by IFisher et al.\ 
(PHH, but less than the value of 32.1° ± 0.7° found bv lStenflo and Kosovichevl 
(|2012p using MDI magnetograms. When we exclude all fields except for 40 kG 
- 50 kG for all the fluxes we consider here, tub increases significantly to 26° ± 
2°, which is closer to the value derived from magnetograms. 

On the other hand, the scatter of the tilt angles around their linear Joy's 
Law fit line is shown to be too large for initial magnetic field strengths of 15 kG 
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and 30 kG with fluxes of 10 21 Mx and 10 22 Mx, as compared to the observed 
value of Fis her et al.\ (|1995p for white light sunspot group images. While the 
scatter of the tilt angle increases with decreasing flux (Table [3]) , we find no 
statistically significant dependence of the Joy's Law trend on flux (Table 0), 
consistent with the results of lStenflo and Kosovichevl (|2012|) . We also find that 
the most common tilt angle produced by our study is 10.6° for tubes with a flux 
of 10 20 Mx, 10 21 Mx, and 10 22 Mx, which agrees with IHowardl (fTOM)) who finds 
that most tilt angles fall within the range of 7.5° - 10° as obtained from Mount 
Wilson magnetogram data, although our average tilt angles are higher. 

Similar to previous studies ( |Caligari et al.\ 119951 IFan and Fisher! 119961 
Web er et aU 120 lTj) . we find that for magnetic field strengths < 50 kG, the leading 
leg of the emerging loop tends to have a larger magnetic field than the following, 
which may provide an explanation for the observed better cohesion of the leading 
polarity of an emerging active region as compared to the following polarity. This 
trend of asymmetry in field strength reverses for tubes with an initial magnetic 
field of > 60 kG. However, it may be the case that the morphological asymmetry 
of sunspot regions is less dependent on magnetic field asymmetry, and is rather 
a result of the retrograde plasma flow inside the flux tube from the leading leg 
into the following leg as suggested by recent simulations of sunspot formation by 
Rcmpel (private communication, 2012). If this is indeed the case, then we can 
not exclude magnetic field strengths of greater than 60 kG from the dynamo- 
gencrated magnetic field regime. A study of the magnetic field of the flux tube 
at the top of the simulation domain suggests typical values of 500 G to 15 kG 
for tubes that reach «21 Mm below the photosphere. 

Observations show that sunspot groups tend to rotate faster than the sur- 
rounding solar surface plasma ( |Howard and Harvey| 119701 
IGolub and Vaianal I1978[) . We use the apparent azimuthal motion of the center 
of the intersections of the emerging loop with a constant r surface near the top of 
the simulation domain as a measure of the rotation rate of the emerging region 
and compare them to the average azimuthal rotation rate of the ASH convection 
simulation at r = 0.95-R©, and the surface rotation rate that we would expect 
assuming the surface shear layer as inferred from helioscismology. For tubes with 
a flux of 10 21 Mx and 10 22 Mx, we find that at high emergence latitudes, the 
average rotation rate of the emerging loops tends to be greater than the inferred 
surface rotation rate for all field strengths considered. At lower latitudes, below 
about 35°, loops with initial field strength >60 kG tend to rotate faster than 
the inferred surface rate, consistent with the observed sunspot rotation rate, 
while loops with initial fields of about 40 kG - 50 kG tend to rotate at a similar 
rate as the surface rate. However, for initial magnetic fields below 40 kG, the 
rotation rate at low latitudes tends to be slower than the surface rate, contrary to 
observations. Thus comparison with the observed sunspot rotation rate seem to 
favor stronger fields as the progenitor of solar active regions. However, because of 
the limitations of our model, we recognize that there remain large uncertainties 
in our results of the azimuthal motion of emerging loops and how they relate to 
the observed sunspot rotations. These limitations arise because our simulations 
stop before the tube enters the solar surface shear layer, and because we do not 
address the deformation and fragmentation of the flux tubes which may result 
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in a stronger and more complex coupling of the magnetic fields with the fluid 
motion, especially in the upper layers of the convection zone. 

Overall, the results in this study suggest that the initial field strength of 
active region progenitor flux tubes needs to be sufficiently large, probably > 
40 kG, in order for them to satisfy the Joy's Law trend for mean tilt angles 
as well as the observed amount of scatter of the tilt angles about the mean 
Joy's Law behavior. Weaker magnetic fields tend to produce too large a scatter 
to be consistent with the observed results. Mid-field strengths of 40 kG-50 
kG, which take the longest to rise, do the best job at matching magnctogram 
observations of Joy's Law dependence. Although only 50 kG or greater field 
strengths can rotate at or faster than the solar surface rate. So, according to 
our thin flux tube approach, magnetic field values need to be of moderate to 
large field strengths for tubes with fluxes of 10 21 Mx and 10 22 Mx to produce 
sunspot rotation behavior. The addition of multiple magnetic flux values to this 
study allows us to perforin a more comprehensive study on the dependence of 
flux tube evolution with regard to magnetic flux. We find, that with convective 
effects, the tube rise time substantially decreases with decreasing flux. There 
also is no statistically significant dependence on the Joy's Law trend in relation 
to the value of the magnetic flux, although the spread of the tilt angles about the 
Joy's Law trend does increase as flux decreases. Convective effects are important 
to flux tube development at these field strengths and magnetic flux values, and 
should be incorporated into future studies. 

With this new study, all magnetic field strengths now show a positive Joy's 
Law trend, and we are now learning that the magnetic field asymmetry of 
sunspots could be due to a retrograde flow of plasma along the emerging flux 
tube. As a result, the estimate of the required > 40 kG magnetic field strength 
is based here mainly on the tilt angle scatter and rotation rates. However, iden- 
tifying the rotation rate of the emerging flux tube, in particular, is an aspect 
where the thin flux tube approach may be lacking in its ability to adequately 
capture all of the physical processes involved. We suggest that this aspect of 
flux tube evolution would be a useful topic for further high-resolution three- 
dimensional MHD simulations of global flux emergence. In addition, more of 
these simulations are needed using realistic active region flux tube parameters 
to further constrain the magnetic field strength generated by the solar dynamo 
based on observational diagnostics. 
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